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Glasses are lucrative engineering materials owing to their superior mechanical properties such as high 
strength and large elastic strain. A central question concerns incipient plasticity - the onset of permanent 
deformation - that is central to their relaxation, aging, yield and fracture. Here, we use an analogue of 
nano -indentation performed on a colloidal glass to obtain direct images of the incipient plasticity, allowing 
us to elucidate the onset of permanent deformation. We visualize the microscopic strain by following 
distorted nearest neighbor configurations, and observe a surprising hierarchical structure of deformation: at 
the onset of irreversible deformation, the strain acquires a robust fractal structure, and we measure its fractal 
dimension. These results give direct evidence that the onset of permanent deformation has the hallmarks of 
a critical point, in agreement with recent theoretical works. 

I ncipient plasticity addresses a crucial step in the deformation of materials. It marks the onset of permanent 
I deformation due to irreversible atomic rearrangements. At small applied strain and short times, solids deform 
I elastically, being able to recover their shape when the strain is removed. However, at larger applied strain and 
on longer time scales, irreversible atomic rearrangements occur, leading to permanent deformation. Glasses are 
frozen on experimental time scales, and the onset of irreversible rearrangements is crucial for their relaxation, 
aging and deformation 1 " 3 . Nanoindentation experiments can provide crucial insight into the incipient plastic 
deformation by resolving individual relaxation events that govern the very earliest stages of plastic deformation 4 " 6 . 
While nanoindentation experiments on conventional materials allow the resolution of individual elementary 
plastic events from force-displacement measurements 7 ' 8 , however, the direct imaging of the elementary relaxa- 
tion, and of the microscopic displacements at the onset of plastic deformation is prohibitively difficult. Colloidal 
glasses provide good models for visualizing the structure and dynamics of glasses in real space. Here, we use an 
analogue of nanoindentation on a colloidal glass to obtain direct images of the incipient deformation. We track 
the individual particles in three dimensions, and directly follow the microscopic strain during the indentation, 
allowing us to elucidate the nature of the incipient deformation. We find robust scaling relations of the micro- 
scopic strain, giving direct evidence that the onset of permanent deformation has the hallmarks of a critical point, 
in agreement with recent theoretical work on fibrous networks 9 . Such critical behavior has important conse- 
quences for the application of amorphous materials as new engineering materials at small scales 1 . 

Results 

Dense colloidal suspensions provide good models to elucidate deformation at small scales; at high density, the 
motion of the particles becomes arrested similar to the atomic motion in conventional glasses. We prepared a 
dense amorphous colloidal sediment roughly 40 particle diameters thick, and used a sewing needle that we push 
slowly into the sediment to indent the colloidal glass (see Methods). The length scale ratio of the needle diameter, 
the colloidal particle size and the film thickness are similar to those used in nano -indentation experiments; hence, 
the experiments are comparable. The needle touches the surface of the sediment at time t=0. The colloidal model 
allows us to follow individual particle displacements deep in the bulk of the colloidal glass. We show particle 
displacements after 21 min of indentation in Fig. lb. Large displacements concentrate below the tip as expected; 
we find that, however, even far away from the tip, pronounced displacements occur. To explore the deformation 
of the glass in more detail, we determine the local strain from distorted nearest neighbor configurations. Nearest 
neighbors of a particle are determined as those separated by less than r l5 the first minimum of the pair correlation 
function. Equally, one could determine nearest neighbors from the Voronoi tessellation of space. We find that the 
result does not depend on the particular choice of nearest neighbors. For each particle, we then follow its nearest 
neighbors in time, and determine the best affine deformation tensor a that transforms the change of nearest 
neighbor vectors over the time interval, At. To do so, the best affine transformation tensor that transforms the 
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Figure 1 | Indentation setup and strain distribution, (a) Schematic showing the needle tip with respect to the amorphous colloidal film, and the imaged 
66 urn by 66 urn by 30 (am volume, (b-d) Particle displacements and strain distribution after 23 min of indentation, (b) Particle displacements in a 5 urn 
thick section centered below the needle. Arrows are enlarged by a factor of two for clarity, (c, d) 66 um by 66 um by 25 um reconstruction of the strain 
distribution. Particle color indicates the magnitude of the maximum local shear strain, e max (see color bar). In (d), particles with e max < 0.1 have been 
omitted for clarity, (e) Maximum shear strain as a function of height below the tip, experimental measurement (red dots) and continuum elasticity 
prediction (dashed curve), (f)-(h) Time series showing the evolution of the maximum shear strain in a 66 (am by 66 (am by 5 (am section, 10 (am below the 
needle, after 9 min (f), 15 min (g), and 21 min (h) of indentation. 



nearest neighbor vectors d^t) at time t, to the vectors d^t+At) at 
time t+At is determined by minimizing the mean square difference 
(di(t+At) — udi(t)) 2 10 . The symmetric part of a corresponds to 
the local strain tensor Sy of the particle under consideration. A good 
measure of the local shear deformation is the maximum shear strain, 
£max = | fii — e 2 1 / an invariant of the strain tensor, where, and e 2 
are the largest and smallest eigenvalues of the strain tensor. The 
maximum shear strain g max reflects the shear strain acting along 
the principal axes of the strain tensor, and has been used as a good 
measure of strain in the inhomogeneous deformation of isotropic 
materials 11 . To visualize its distribution, for each particle, we deter- 
mine values of g max from the eigenvalues of the strain tensor assoc- 
iated with the particle. Reconstructions of the resulting strain 
distributions are shown in Fig. lc and d. Orange and yellow spheres 
indicate particles with large £ max , while blue particles indicate small 
e max . In Fig. Id, we have omitted particles with e max <0.1 1 for clarity. 
Orange particles accumulate below the tip, indicating strong shear 
strain. A strain maximum occurs at roughly 10(am below the tip 
(Fig. le) in good agreement with continuum theory predictions. In 
contrast to continuum elasticity predictions, however, large shear 
strain occurs surprisingly far away from the tip. To demonstrate this, 
we focus on a horizontal section, 10 (am below the tip, and show the 
time evolution of strain in Figs. If to h (see supplementary video). 
Already at early times (Fig. If and g), large shear strain occurs 
throughout the field of view as indicated by the occurrence of yellow 
spheres far away from the tip. This leads to an unexpected distri- 
bution of strain at later times, as demonstrated in Fig. 2a. Large shear 
strain concentrates at preferred distances to the origin (dashed cir- 
cles), indicating a cascade-like structure. Similar cascade structures 
are observed in all components of strain. As an example, we show the 
distribution of normal strain s zz in Fig. 2b. While blue spheres in the 
center reflect the compression induced by the tip, the pattern of blue 
and red spheres indicates a cascade of compression and dilation 
similar to that of the maximum shear strain (Fig. 2a). Obviously, this 
cascade structure is a characteristic property of the emergent plastic 
deformation; to show this, we compare the observed strain 
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Figure 2 | Cascade structure of the incipient deformation, (a, b) Strain 
distribution in a 5 (am thick section 10 (am below the tip after 21 min of 
indentation, a, Distribution of the maximum shear strain, e max . 
b, Distribution of normal strain e zz . (c) Angle- averaged strain as a function 
of distance from the center of the tip. Symbols indicate the measured values 
of £ ma x (red dots), e zz (green dots), and the radial shear strain e rz (violet 
dots). Dashed lines indicate continuum predictions for e max > e zz> an d £ rz > 
respectively. Arrows indicate pronounced deviations from the predicted 
elastic distribution associated with the onset of permanent deformation. 
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Figure 3 | Time evolution of strain, (a)-(c) Time series showing the evolution of regions of high strain after 16, 20, and 24 min of indentation. Red 
arrows indicate highly deformed zones penetrating the sample with time, (d) Time correlation of strain between the inner and middle (blue), middle and 
outer (red) and middle and outer (green) dashed circle in Figs. 3a-c. The time correlation is calculated by correlating the average strain within cylindrical 
shells of thickness 4 (am, and radii Ri= 10, R 2 =25 y and £3 = 38 (am, centered at the origin. Maxima indicate the characteristic delay times of plasticity, 
Ati2 = 6 y At 2 3 = 8 and At i3 =14 min. 



distribution with the elastic strain distribution predicted by con- 
tinuum elasticity. We assume the glass to be an isotropic elastic 
medium with Young's modulus, E, shear modulus, fi, and Poisson 
ratio v=0.34 12 . The elastic response of an indented elastic half- space 
has been treated analytically 1314 ; the elastic strain depends only on 
the radius of contact, ,R C , between the indenter and the medium, and 
the maximum normal strain s 0 in the center of contact. We measure 
R c directly from the imaged contact surface between the tip and the 
colloidal glass and obtain R c =22 (im. We then adjust £ 0 m order to 
best fit the experimental strain distributions, and obtain 8 0 = 0.055. 
The resultant elastic strain distributions are shown in Fig. 2c (dashed 
lines). Also shown by dots are the measured strain distributions. 
Striking differences occur at Ri = 10, R 2 =25 and #3 = 38 |im (arrows) 
in all components of strain. These distances coincide precisely with 
the location of the circles in Fig. 2a and b; we associate these devia- 
tions from the elastic strain distribution with the onset of plastic 
relaxation. 

Load- displacement curves of metallic glasses exhibit characteristic 
bursts at the onset of plastic flow 4 that have been associated with the 
nucleation of shear bands, where slip occurs collectively along a 
surface. Such shear bands occur at an intermediate regime of 
indentation rates where the glass exhibits inhomogeneous flow 8 . 
The long time scales afforded by our colloidal glass permit us to 
access the regime of uniform flow at low indentation rates 15 : 
Surprisingly, we observe a robust cascade structure of the incipient 
plastic deformation, demonstrating the occurrence of strong correla- 
tion in this uniform flow regime. 



The long time scale of our colloidal glass allows us to even follow 
the evolution of this cascade in time: We focus on regions with large 
shear strain, and follow their time evolution in Fig. 3a-c. The incipi- 
ent plastic deformation penetrates the glass with a characteristic 
propagation time. To determine this time scale, we measure the time 
delay between the plastic activity at the distances R i9 R 2 and £3 in 
Fig. 2. We correlate the angle-averaged strain e,- (t) at R t at time t with 
the strain 8j(t + At) at Rj at some later time, t + At. The correspond- 
ing correlation functions, shown in Fig. 3d, reveal distinct time cor- 
relation: Maxima at At=6, 8 and 14 min demarcate characteristic 
delay times associated with the propagation of plastic deformation, 
and reflect the time necessary to form new plastic zones in the strain 
field of existing ones. Thus, along with the spatial correlation, there is 
also strong time correlation in the incipient plastic deformation. 

This strong spatial and temporal correlation results in a complex, 
highly correlated network of deformation. To elucidate it, we explore 
the analogy with fractal structures encountered in the aggregation of 
particles, where strong interparticle attraction leads to hierarchical 
structures characterized by a robust fractal dimension 16 . While for 
space-filling structures, the number N(r) of particles within a dis- 
tance r grows as N(r) oc r°, where D is the dimension of space, in 
open, fractal structures, N(r) oc r df , where dj<d is the fractal dimen- 
sion. Surprisingly, we find that a similar description applies to the 
distribution of strain during the incipient deformation. To show this, 
we again focus on the high-strain particles shown in Fig. 3a-c, and 
determine their number N(r) within a lateral distance r to the origin. 
We plot AT as a function of r in Fig. 4a, inset. Within the highly 




Figure 4 | Fractal deformation field, (a) Scaling of the number N of high-strain particles, normalized by the homogeneous distribution, as a function of 
distance r from the origin. Two regimes with different scaling are delineated by the vertical dashed line. For r< 15 (am, N is constant, indicating a compact, 
space-filling structure, while for r>15 [imN ~r~ 12 indicating a fractal structure. Colors indicate increasing time intervals of indentation, t=l5 (red), 16 
(violet), 17 (gray), 18 (green), 19 (blue), 20 (yellow) and 24 min (turquoise). Inset: Same data without normalization shows the fractal scaling of the 
number of high- strain particles according to N(r) ~ r 2 , and N(r) ~r 0-8 . (b) Relative frequency of strain values as a function of strain magnitude for 
particles with distance r< 10 (am (blue) and r>20 (am from the center (red). The dashed lines indicate fits to a Gaussian (red) and a power-law distribution 
(blue). 
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deformed center, iV—r 2 , indicating a compact space-filling structure. 
Outside the central zone, however, N grows slower, with a slope of 
— 1.2 indicating a fractal structure of deformation. To show this most 
clearly, we normalize by the homogeneous distribution, N = npr 2 , 
where p is the average density of high-strain particles, and show 
the normalized number of particles in Fig. 4a, main panel. The 
change from a horizontal to a negative slope of the data indicates 
the transition from a compact to a fractal structure of deforma- 
tion. The excellent collapse of the data taken after different times 
of indentation demonstrates that this fractal structure is a robust 
property of the incipient deformation of the glass. This fractal 
structure of deformation allows plasticity to penetrate the glass 
at a low density. While our field of view limits the range of acces- 
sible length scales, we obtain further evidence of the fractal nature 
of deformation from the probability distribution of strain. We 
show histograms of strain values in Fig. 4b, where blue and red 
dots indicate particles in and outside the central zone, respectively. 
The data reveals the very different character of the distributions: while 
inside the central zone, particles exhibit a Gaussian strain distribution 
(blue dots and dashed line), indicating plasticity occurs in a statistically 
independent manner, outside the central zone, particles exhibit a 
power-law distribution (red dots and dashed line), indicating strong 
correlation, in agreement with the fractal structure of deformation. We 
hypothesize that due to the high shear stress in the center, plastic zones 
can form independently of each other, while at smaller stresses further 
away, plastic zones can only form when triggered by adjacent plastic 
zones, and hence plasticity occurs in a highly correlated fashion. These 
strong correlations result in a hierarchical structure of deformation. 

Discussion 

Recent theoretical work on disordered fiber networks 9 indicates that 
the loss of network rigidity is associated with a true critical point. 
Our observation of a hierarchical strain distribution indicates that 
even a glass shows hallmarks of critical behavior at the loss of 
elasticity and the onset of plastic flow. Critical points demarcate 
qualitative changes of behavior, typically associated with symmetry 
changes. For the incipient deformation investigated here, the sym- 
metry change occurs in time, and demarcates the transition from 
reversible elastic to irreversible plastic deformation. While the direct 
imaging of microscopic strain allows us to visualize its fractal struc- 
ture in a colloidal glass, we expect our results to be general and apply 
to molecular glasses as well. We estimate from approximate con- 
version of the deformation rate from indentation to uniaxial strain 
rate that the strain rate of our experiment is ~3X 10~ 5 s~\ a factor 
of three larger than the inverse structural relaxation time of our 
glass. This ratio is typically encountered in the high-temperature 
deformation of molecular glasses. For example, for typical strain 
rates employed in nanoindentation experiments of atomic glasses 
in the range of 10" 3 s" 1 to 1 s~\ the same ratio of imposed strain 
rate to relaxation rate is found in the range between 0.9 and 1.1 of 
the glass transition temperature. Deformation tests on metallic 
glasses show that indeed homogeneous deformation occurs at these 
strain rates and temperatures 8 , in agreement with our experiments. 
Therefore, our observations should be relevant for these atomic 
glasses as well. 

The application of indentation on a colloidal glass reveals a sur- 
prising hierarchical structure of strain at the onset of permanent 
deformation. This hierarchical structure is reminiscent of a sec- 
ond-order transition, indicating that the glass at the onset of flow 
exhibits hallmarks of critical behavior, in agreement with recent 
theoretical work on fiber networks. The similarity between the fractal 
structure of deformation and that encountered in fractal aggregates 
of particles suggests an intriguing analogy between the mechanisms 
of structural arrest and shear-induced relaxation. Finally, the strong 
correlation of the incipient strain field should have important con- 



sequences for the deformation in small dimensions, such as the one 
exploited in applications of amorphous materials at small scales 1 . 

Methods 

We prepare a colloidal glass by sedimenting 1.55 um diameter silica particles onto 
a cover slip by centrifugation. The particles are suspended in a mixture of water 
and dimethylsulfoxide (DMSO) that matches their refractive index. The cover slip 
surface is coated with a 5um thick layer of polydisperse colloidal particles to avoid 
boundary-induced crystallization. The resulting colloidal amorphous film has a 
volume fraction of <^~0.60, larger than that of the colloidal glass transition, 
0 g ~O.58 17 , and is therefore well inside the glassy state. A sewing needle attached to 
a piezoelectric translation stage is used to indent the colloidal glass. The needle 
has an almost hemispherical tip with radius —38 urn, and produces a strongly 
inhomogeneous stress field with a distinct maximum well below the contact 
surface. We lower the needle slowly at a rate of 2.9 um/h into the amorphous 
sediment and use confocal microscopy to image individual particles in a 66 um by 
66 um by 48 um section of the glass below the tip (Fig. la). This technique allows 
us to follow the incipient plastic deformation directly at the particle scale. We 
determine individual particle positions in three dimensions with an accuracy of 
0.03 um in the horizontal, and 0.05 um in the vertical direction 18 . We acquire 
image stacks every 60 sec over a time interval of 25 min to track the motion of 
—3X10 4 particles in the imaged section. One image needs 45 sec to acquire. 
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